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Abstract. In order to investigate photoinduced insulator-metal transitions observed in 
correlated electron systems, we propose a new theoretical method, where we combine a Floquet- 
matrix method for AC-driven systems with the dynamical mean-field theory. The method 
can treat nonequilibrium steady states exactly beyond the linear-response regime. We have 
applied the method to the Falicov-Kimball model coupled to AC electric fields, and numerically 
obtained the spectral function, the nonequilibrium distribution function and the current- voltage 
>-^J , characteristic. The results show that intense AC fields indeed drive Mott-like insulating states 

C"| ' into photoinduced metallic states in a nonlinear way. 

O ' 

J> ' 1. Introduction 

There is an increasing fascination with photoinduced phase transitions (PIPTs), where a 'phase 
transition' is triggered by a laser of strong intensity. Specifically, strongly correlated electron 
systems harbor the phenomena, where an insulating state in, e.g., manganese oxides can be 
made metallic when a pulsed laser is injected [U In some cases a transient ferromagnetic 
spin alignment can be generated as observed in recent pump-probe spectroscopy experiments [3j. 
While these phenomena are intrinsically nonequilibrium, a remarkable feature is the threshold 
behavior, namely, the transition only occurs when the intensity of light exceeds a certain 
threshold value. This implies that PIPT is essentially beyond the linear-response regime for 
the external driving fields. It is then a theoretical challenge to identify the nature of the states 
emerging in nonequilibrium against the frequency and the intensity E of the AC field. 

A crucial aspect in theoretically studying PIPT in strongly correlated systems is that we 
have to simultaneously take account of both the electron correlation and the nonlinear effect of 
the electric field. Here we propose a theoretical method, in which we incorporate the Floquet- 
matrix method for AC field problems [H [5] into the nonequilibrium dynamical mean-field theory 
(DMFT) [6]. The Floquet approach can treat a nonequilibrium steady state nonperturbatively 
{i.e. up to arbitrary orders of E), while the DMFT can describe correlation effects like Mott's 
transition. Here we first sketch the method in section [21 and then apply it to the Falicov- 
Kimball model (one of the simplest model of correlated electrons that can be solved exactly 
within DMFT) in section [3l From the obtained results, we discuss how the Mott-like insulator- 
to-metal transition is induced by AC fields in section [H 
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2. Floquet-matrix method combined with DMFT 

We outline how we can incorporate the Floquet-matrix method in the DMFT. Details of the 
method are described elsewhere [3 [8]. Our argument is based on the Keldysh Green's function 
formalism for nonequilibrium states. In general, a Green's function, G(t,t'), for a system out of 
equilibrium should depend on two time arguments t and t' independently (while in equilibrium 
it depends only on t — t'). The original idea of the Floquet method is that, when the driving 
field is periodic in time with the period r = 2ir/Q, there is a time-analog of Bloch's theorem 
for spatially periodic potentials. One consequence is that G satisfies a periodicity condition 
G(t + T,t' + r) = G(t, t'). Then one can perform the Wigner transformation on G, 
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dt rcl - dt avc e^ +inm -G(t,t'), (1) 
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where t re \ = t — t' and t ave = (t + 1')/2. With G n (uj) we define the Floquet-matrix form of G as 



u H - — Q J . (2) 

For G mn (uj) we use the reduced zone scheme, i.e., the range of uj is restricted to the 'Brillouin 
zone' on the frequency axis, —Q/2 < to < £7/2. With this setting G n {uo) has a one- 
to-one correspondence to G mn {uj). An advantage of employing the form (J2j) is that the 
convolution of two functions J dt'A(t,t')B(t',t") is mapped to the product of two matrices 

Yin' A n n'(u)B n > n ii (u). 

We can then employ the definition ([2]) to rewrite all the self-consistent equations in DMFT 
in a Floquet-matrix form, which can be solved by the numerical iteration procedure as in the 
equilibrium case. The fact that we can take advantage of taking the inverse of a matrix when 
solving the Dyson equation makes our calculation numerically quite efficient. 

As an input for DMFT we need the inverse of the retarded Green's function for noninteracting 
electrons, which is given by 

( G k° )mn(u) = (uj + nU + /J, + ir])S mn - (e h ) mn , (3) 

where [i the chemical potential, rj = + , and A(t) the vector potential. We assume that eA(t) = 
eE sin(f2t) / Vl [1, 1, . . . ,1] (representing a homogeneous AC electric field in d dimension). In the 
hypercubic lattice, (efc) mn = ekJ m - n (eE/Q) (m — n: even), or itkJ m -. n [eE /Q) (m — n: odd), 
where J n is the Bessel function, e& = — 1/VdYli cos ^ an d £fc = — l/v^X^ s ™ ^ n integral 
over k is calculated with the joint density of states p(e, e) = e~ e ~ e /it [6]. 



3. Falicov-Kimball model driven by AC electric fields 

Having constructed the theoretical framework, we are now in position to demonstrate its 
application to the Falicov-Kimball model with a heat bath. The Hamiltonian is -fftot = 

H sys + -ffsys-bath + -^bath w ith 

= e k-eA{t)C ] k C k + U J2 C i C ififi> ( 4 ) 
fc i 

where c\ creates an itinerant, spinless electron while fj creates a localized one, -ffbath represents 
the bath, and -ff S ys-bath the coupling between the system and the bath. We take a noninteracting 
fermion bath attached to each lattice site in the system. After the bath degrees of freedom are 
integrated out, the correction to the self-energy T(uj) is a constant, T, independent of uj in the 
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Figure 1. The local spectral function Aq(uj) (solid curves) and the particle number density 
Nq(lo) (shaded regions) (left panels), and the effective distribution function / e ff(w) (right) for 
E = (a), E = 1.0 (b), E = 2.0 (c) with VL = 1.0, U = 3.0, T = 0.05 and T = 0.05. 




Figure 2. The gauge invariant spectral function AQ{k,u)) (color-coded), and the local spectral 
function Aq{uj) (curves) for E = (a) and E = 0.6 (b) with Vt = 1.0 and U = 2.2 on the three 
dimensional cubic lattice, and k = k(l, 1, 1). 



wide band limit. Assuming that the bath is in equilibrium with a temperature T, we can define 
the nonequilibrium steady state (parametrized uniquely by T and T). Note that the steady 
state is determined without the noninteracting Keldysh Green's function G K0 . Namely, the 
information on the initial state of the system and the way to switch on the field, on which G K0 
depends, is wiped out due to the dissipation caused by the bath. 

We consider the following quantities: the spectral function A n (uj) = — ^ fc Im[(G^) n (a;)]/7r 
(FigHJ left panels), the particle number density N n (co) = — i Ylki^k )n( tJ )/2 7r (FigH])) the 
effective distribution function f c ^{uj) = Nq(u)/Aq (uj) (FigHJ right), the gauge invariant spectral 
function A n (k,u;) = -Im[(G%) n (uj)]/ir (FigOJ), and the current j(t) = ie^2 k e k _ eA ^G^(t,t) 
(Figj3]). Here G R and G < are the retarded and the lesser Green function respectively. Note 
that, to make the fc-dependent Green's function gauge invariant, we use the modified function 
G(t,r;t',r') = exp(ief^' ) dx»A u (x))G(t,r;t',r') [9]. 
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Figure 3. The current-voltage characteristic (solid curve) with f2 = 1.0, U = 3.0, V 
T = 0.05. We also show the linear-response result (dashed line) for comparison. 
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4. Discussion 

From the results one can see how the Mott insulating state is driven into a metallic state as one 
increases the intensity of the AC field beyond the linear-response regime. Specifically, we find 
in Fig(T]that the Mott gap collapses, and a new spectral weight emerges in the midgap region. 
Since the optical gap ~ U(= 3.0 here) is greater than the frequency f2(= 1.0 here), the direct 
interband transition is forbidden. However, if the intensity E is sufficiently large, the electrons 
are excited to the conduction band through the midgap state, and the distribution function 
/eff(w) (Fig[T]) considerably deviates from the Fermi distribution f(uj) = l/(e w l T + 1) realized in 
equilibrium. More surprisingly, there is a region around certain uj where the population exceed 
those in lower-energy regions, a kind of population inversion. In the band dispersion of the 
system out of equilibrium shown in Fig|2l we again observe that the midgap state is induced by 
the field, whose band structure is roughly replicas of the original one with the spacing £1. Finally 
in Fig|3j the nonlinearity with the AC field is captured most clearly as the I-V characteristics 
obtained in the present method. Here J = Re - f^dt e int j(t) is the current. When E is small, 
the current is approximately proportional to E as in the linear-response theory, but a nonlinear 
response becomes evident for E > 0.2. This suggests that in the heart of the photoinduced 
insulator-metal transition lies the nonlinearity of the driving field. 

During the preparation of this paper, we notice that a similar idea of using a matrix form 
was mentioned in ref. |10j . 
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